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Estimation  and  Control  of  Parameters 
in  Linear  and  Nonlinear  Distributed 
Models  of  Flexible  Structures 

Final  Report 
February,  19S4 
1.  Introduction  and  Goals. 

This  project  seeks  to  study  the  estimation  of  elastic,  damp¬ 
ing,  and  material  parameters  in  flexible  structures.  Of  particu¬ 
lar  interest  are  problems  in  the  design  and  estimation  of  param¬ 
eters  in  structures  made  up  of  systems  of  coupled  beams  and 
plates,  the  estimation  of  parameters  in  models  that  may  not  have 
unique  solutions,  and  the  estimation  and  design  of  various  plate 
and  shell  models  incorporat ing,  for  example,  large  deformations, 
variable  thickness,  existing  curvatures,  contact  and  possibly 
friction  conditions. 

Because  of  their  application  as  constructional  elements,  it 
seems  important  to  study  both  the  static  and  dynamic  behavior  of 
beams,  plates,  and  shells,  as  well  as  their  combinations.  Accord¬ 
ingly,  this  project  seeks  to  study  the  identification  of  param¬ 
eters  and  the  design  of  structures  modeled  by  distributed  systems. 
Such  efforts  are  fundamental  to  our  ability  to  simulate  and  con¬ 
trol  their  behavior.  Our  goals  are  (i)  develop  and  test  numerical 
algorithms  for  the  estimation  of  parameters  and  design  of  certain 
structures  composed  of  connected  beams  and  plates,  <ii)  develop 
theories  and  test  numerical  algorithms  for  the  estimation  of 
parameters  in  models  that  may  have  multiple  states,  (iii)  use  the 
above  theories  to  consider  various  models  of  plates,  shells,  and 
beams  with  possible  nonlinearities,  and  in  cases  where  feasible 
(iv)  test  algorithms  with  laboratory  data  when  possible. 

An  identification  or  estimation  problem  seeks  to  determine 
parameters  within  a  mathematical  model  from  observed  data.  The 
central  issue  is  how  to  utilize  the  data  to  determine  the  desired 
parameters  within  the  context  of  the  model.  A  design  problem 
views  parameters  within  the  problem  as  controls  to  be  adjusted  so 
as  to  produce  desired  results  under  certain  conditions.  Although 
these  problems  are  stated  differently,  a  common  mathematical  opt¬ 
imization  formulation  may  be  given  to  them.  Available  is  a 
model  equation 
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<1.1)  C<q)  u<q)  *  f 

in  which  the  parameters  q  to  be  estimated  belong  to  a  specified 
admissible  set  Qad  of  a  Banach  space  Q  and  solutions  u  *  u(q) 
belong  to  a  Banach  space  X  of  states.  The  set  Qad  should  be 
physically  meaningful.  However ,  whereas  in  previous  work  it  has 
been  a  condition  that  the  parameter-to-state  mapping  is  well- 
defined  as  a  function  from  Qad  into  the  state  space  X  we  will  in 
this  work  allow  the  possibility  of  multiple  solutions  for  the 
st  at  e-equat i on . 

Available  (for  the  estimation  problem)  are  data  z  that  we 
view  as  belonging  to  an  observation  space  Z.  An  approach  to 
solving  the  estimation  problem  is  formulated  as  a  minimization 
problem! 

Find  q0«  Qa:?  such  that 

(1.2)  J(qQf z)  -  infCJCq, z)!  q  «Qad> 

admissible  functions.  An  example,  if  sense  can  be  given  to  it,  of 
a  fit-to-data  functional  J(q,z)  is  the  so  called  regularized 
output  least  squares  functional 

<1.3)  J(q,  z )  ■  ||Cu(q)-z||£  +  p  llq  11^ 

where  Z  is  the  observation  space  (Hilbert  space),  C  is  the  obser¬ 
vation  operator  that  takes  X  into  Z,  and  p  >  0  is  a  regulariza¬ 
tion  parameter.  If  there  may  be  multiple  solutions  to  (1.1),  care 
must  be  given  in  the  formulation  of  the  fit-to-data  functional 
such  that  it  is  a  well-defined  functional  that  is  bounded  below, 
is  weakly-lower  semicont inuous  over  Qad  with  respect  to  a  topology 
in  which  Qad  has  suitable  compactness  properties.  There  are  other 
approaches  to  identification  but  these  typically  assume  the 
parameter-to-state  mapping  is  well-defined  and  we  refer  the  reader 
to  Cl]  for  discussions.  Issues  of  interest  for  (1.1)-(1.2) 
include  the  proper  formulation,  regularity  properties  of  solutions 
of  (1.2),  approximation,  stability  with  respect  to  the  data,  and 
uniqueness  or  identif iabil ity. 

The  minimization  problem  (1.2)  may  also  be  viewed  as  design 
problem  where  z  represents  a  target  function  that  we  wish  to  reach 
by  choosing  a  suitable  parameter  q.  Certainly,  it  makes  sense  from 
a  design  point  of  view  to  impose  as  a  condition  of  the  design 
problem  that  the  mapping  q  *->  u(q)  is  well-defined,  although  one 
does  not  have  to.  Design  problems  seem  to  be  less  data  oriented 
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than  estimation  problems  sines  on*  can  of tan  exercise  more  control 
over  the  formulation.  In  estimation  problems  one  is  limited  by 
the  ability  to  measure  the  state.  The  data  are  not  always  in  the 
most  desirable  form.  Indeed  this  is  one  of  the  first  lessons  that 
one  learns  attempting  to  use  experimental  data. 

If  one  expects  to  use  estimation  techniques  to  identify 
parameters  from  measurements,  then  one  must  first  decide  upon  the 
a  mathematical  model  (1.1).  Hence,  one  must  decide  whether  the 
mathematical  model  is  appropriate  for  the  physical  system  being 
observed.  The  mathematical  model  embodies  the  pertinent  princi¬ 
ples  and  assumptions  based  on  physics  and  continuum  mechanics. 
Certainly,  one  must  be  cautious  about  using  models  to  fit  data 
that  do  not  satisfy  the  assumptions  upon  which  those  models  are 
derived.  It  is  always  desirable  to  validate  in  some  way  the  model 
in  question.  For  example,  comparing  spectral  data  predicted  by 
the  model  with  spectral  data  observed  in  the  laboratory  provides  a 
reasonable  step  in  model  comparison. 

Having  proposed  the  mathematical  model,  one  must  address  how 
to  use  the  data  within  the  framework  it  imposes  for  the  purpose  of 
estimation  and  mode  validation.  The  data  for  systems  of  interest 
in  this  project  are  obtained  as  "pointwise"  measurements  of 
deformation  or  strain  in  the  static  case  or  acceleration  or  vel¬ 
ocity  or  their  Fourier  transforms  at  points  in  the  dynamic  case. 
These  data  should  be  in  a  form  consistent  with  the  mathematical 
model  having  the  proper  units,  etc.  The  mathematical  analysis  of 
identification  problems  such  as  (1.1>-(1.2>  gives  information  on 
how  best  to  formulate  estimation  problems  to  use  these  data. 

In  the  following  we  indicate  some  of  the  results  we  have 
obtained  during  the  last  year.  In  Section  2  we  discuss  the 
estimation  of  elastic  parameters  in  dynamic  von  Karman  plate 
models.  In  Section  3  we  report  on  our  work  with  Prof.  D.  L. 
Russell  on  the  modeling  and  spectral  identification  for  so  called 
narrow  plate  models  that  should  be  useful  in  modeling  wings,  fan 
blades,  and  corrigated  structures.  Finally,  in  Section  4  we 
indicate  our  work  on  the  modeling  and  estimation  of  parameters  in 
dynamic  models  of  coupled-beam  systems. 

2.  Estimation  of  elastic  parameters  in  a  dynamic  nonlinear  plate 
model . 

During  the  funded  period  another  ther  area  of  investigation 
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has  bean  the  estimation  of  parameters  in  static  and  dynamic 
versions  of  plate  models  that  allow  large  deformations.  For  the 
dynamic  problems  we  have  considered  the  von  Karman  and  Morozov 
models.  It  is  well-known  that  global  solutions  of  the  von  Karman 
equations  are,  in  general,  nonunique  £23  although  local  uniqueness 
of  classical  solutions  for  certain  cases  has  been  demonstrated  in 
£3,43.  In  order  to  give  a  mathematical  formulation  of  the  param¬ 
eter  ident i f ication  problem  for  either  the  static  or  dynamic  von 
Karman  equations,  it  is  necessary  to  study  the  estimation  of 
parameters  for  systems  that  may  not  have  well-defined  parameter- 
to-state  mappings.  We  thus  consider  a  theory  for  the  estimation 
of  parameters  for  set-valued  par ameter-to-state  mappings. 

We  distinguish  between  two  cases  static  and  dynamic.  For  the 
static  problem  it  is  apparent  the  the  Galerkin  systems  may  not 
have  unique  solutions.  Our  approach  is  to  introduce  model  error 
equations.  See  £83  and  the  description  of  the  model  error  method 
in  our  original  proposal .  One  advantage  that  we  have  discovered 
with  the  model  error  method  in  the  static  case  is  that  we  may  give 
formulations  to  avoid  having  to  calculate  the  adjoint  completely 
and  may  be  able  to  compute  the  derivatives  of  the  fit-to-data 
functional  directly  from  the  state  solutions.  This  obviously  can 
lead  to  significant  computational  savings.  Alternatively,  for  the 
static  case  we  may  seek  conditions  that  imply  that  the 
par ameter-to-state  mapping  is  well-defined.  In  £93  this  approach 
is  taken  that  is  based  on  the  contraction  mapping  principle  for 
existence  as  opposed  to  the  Schauder  fixed  point  principle  £23. 

In  the  dynamic  case  we  can  take  advantage  of  fact  that  the 
Galerkin  system  forms  a  system  of  ordinary  differential  equations 
for  which  we  may  obtain  the  existence  of  unique  solutions.  Our 
approach  is  first  to  formulate  and  analyze  the  parameter  identi¬ 
fication  problem  for  the  von  Karman  problem  using  an  existence 
theory  for  set-valued  par ameter-to-state  mappings.  A  formulation 
is  given  that  utilizes  the  weak*  compactness  of  the  solution  set 
for  admissible  parameters.  There  are  complicating  technical 
factors,  however,  that  imply  different  convergence  results 
£10,113.  We  show  existence  of  a  solution  to  the  output-least- 
squares  estimation  problem.  Further,  we  obtain  certain  weak 
stability  results  and  convergence  properties  for  the  estimation 
problems  of  associated  Galerkin  systems.  Next  we  use  the  fact  that 
the  Galerkin  system  possesses  a  par ameter-to-state  mapping  that 
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is  well-defined  to  pose  well-defined  out put -least -squares  estima¬ 
tion  problems.  We  then  determine  regularity  and  approximation 
results  for  the  optimal  estimators  of  the  Galerkin  problems  and 
use  them  to  prove  stability  and  convergence  results  for  the 
estimation  problems  for  which  the  Galerkin  systems  are  the 
underlying  systems. 

We  give  the  estimation  problem  for  the  Galerkin  problem  and 
report  the  results  of  several  numerical  experiments.  The  results 
of  other  numerical  studies  as  well  as  the  theoretical  development 
are  reported  in  references  C103  and  C 1 1 3 .  We  recall  the  von 
Karman  system 


(2.1)(i)  Uj.£  +  a  <*  a  u)  +  e  CB(u,u),u3  *  f  in  Q 
with  initial  conditions 

<2. 1 ) < i  i  ) 

ut<0)  *  uj 

and  boundary  conditions 

<2.1)<iii>  u  =  |}f=Oon  <0,T> 

where 

<2. 1) <iv) 


a?  B(u,v)  ■  Cu,v3 


with  clamped  boundary  conditions.  Let  k=»l,...,N>  be  linearly 

independent  subset  of  functions  in  V  and  {yk:  k=l,...,M>  be  a 

linearly  independent  subset  of  functions  in  Q.  We  denote  by  VN  = 
spanC^s  k-l,...,N>  and  QM  *  spant^i  k=l,...,M>.  We  set  uN  = 

N  m  M 

,Z.ck<t)  and  an  =  .£  a*  y*.  Defining  the  matrices 
k*l  K  K  i*>l  11 

<Q0>ij  =  4  *iTf  dx  and  <G2>U  "  SQ  *  * i  *  Tj  d* 

<Qk)iJ  "  JQ  \  A  d*  for 

and 

<Hk)ij  *  JQ  EVi»  \  d*  for  k-lf...,N, 

we  set 


H<c,d) 


d 


H * 


Wd)  c 


Finally,  we  define  the  N  x  N  matrix 
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M  u 

GCa)  ■  Z  Gk 

k»l  k 

The  following  semi-discrete  system  is  obtained 

60  ctt  +  GCa)  c  +  c  HCc ,  d)  -  f 

<2. 2)  G2  d  *  HCc,c> 

G0  c<0)  * 

G0  ct<0)  - 

We  use  the  following  difference  equations  for  a  time  approxi- 
mat ion 


i  +  1 


2  c. 


*  ci-l 


+  GCa) 


■i+1 


+  2  c. 


ci-l 


+  c  HCcifd|)  «  fj 

where  h  denotes  the  time  difference  and  the  subscript  i  indicates 
the  time  ti  ■  i  *  h.  Thus,  the  following  sequence  is  generated 

Q0  c0  "  ->*0 

Q0  C1  "  *0  +  h  *1 
G2  dj  -  H<c j ,  ct) 

2 

<G0  +  C§02  GCa))  ci  +  1  »  2  (  G0  -  J  GCa))  c*  - 
-  CG0  +  C§-)2  GCa))  cj.j  +  h2  C  ft~  c  HCcifdt)) 


The  derivative  system  is  obtained  by  formally  differen¬ 
tiating  the  equations  in  C2.1)  with  respect  to  a q  for  8  = 

1, . . . ,M.  Thus,  we  see  that  the  partial  derivative  of  c  and  with 
respect  to  a^  D  jc  and  Dgd,  satisfy  the  system 

G0  CDpO”  +  CGCa)  +  c  Wd))  CDgO  + 

♦  c  WO  D  gd  -  -  Qf(S  c'  +  c) 

<2*3)  G2  CDgd)  -  2  Wc)  CDgO 

with  initial  conditions  given  by 

CD  gO  CO)  -  CDgc)f  CO)  -  0 

To  specify  the  fit-to-data  functional,  we  define  the  MxM 
matrix  G*  by 

CG^ij  -  <*if  Vv 

for  i,j  ■  the  N-vector  valued  function  t  *— >  eCt)  by 
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«i<t)  »  C  »  z)H 

for  i«l,...,N,  and  the  real  valued  function 

c2<t>  -  JQz2(x,t)  dx. 

With  these  definitions  the  fit-to-data  functional  becomes 
<2. 4)  J  (a)  *  jJ(c*(t)  G0  c  ( t )  -  2  e*(t)  c(t>  +  cz(t>>  dt  ♦ 

+  p  a*  G*  a. 

Thus,  J  is  a  mapping  from  rf1  into  R.  We  may  calculate  the  Frechet 
partial  derivative  with  respect  to  a p  to  obtain 

DpJCa)  «  2  fo  (  G0  c(t)  -  «(t))*  (DjcXt)  dt  +  2  P  (a*  G*>  Q 

where  D  pc  is  the  solution  of  (2.3). 

To  solve  the  estimation  problem,  we  attempt  to  minimize  the 
functional  J.  Our  approach  is  to  use  a  steepest  descent  method  to 
obtain  a-vectors  that  decrease  the  size  of  J(*).  We  present  the 
results  of  several  numerical  experiments.  Our  computations  are 
conducted  for  Q»  (0,1)  x  (0,1)  on  which  we  specify  a  mesh  ob¬ 
tained  from  uniform  meshes  on  (0,1)  Aj  with  8  subintervals  and 

a 2  with  3  subintervals.  The  basis  functions  ■C  ^  i*««i  arB  obtained 
as  tensor  products  of  cubic  B-spl ines  over  the  mesh  ^  that  are 
adjusted  to  satisfy  the  clamped  essential  boundary . condit ions. 
Hence,  for  our  examples,  N  =  49.  On  the  other  hand  for  the 
approximat ion  of  the  parameter,  we  use  cubic  B-spl ine  over  the 
mesh  *2  with  no  restriction  on  the  boundary.  The  basis  functions 

(^i>^«i  are  obtained  as  tensor  products  of  these  functions.  Hence, 
M  «  36. 

We  consider  the  following  test  problems.  Let  c  *  0.1.  Defin¬ 
ing  the  function 

uT(x)  -  16  x2  (1  -  x)2, 
we  specify  the  deformations  as 

u^s^(x,y,t)  *  uj(x)*uj(y)*cos(t)  in  S3  x  (0,T). 

Hence,  the  initial  conditions  for  our  examples  are  given  by 

u^(x,y)  ■  uT(x)  *  uT(y)  in  Q 

and 

uj  (x,y)  *  0  in  8. 
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To  carryout  various  experiments  to  recover  the  coefficient  a,  we 
specify  a  particular  coefficient  *otst*  w*  generate  the 

resulting  force  vector  by  using  the  weak  form  of  the  equations 
(2.1)  with  basis  functions  9j  replacing  <P.  Thus,  from  the  data 
u^t  as  an  observation  and  f  as  an  assumed  known  force,  we  wish  to 

recover  the  coefficient  a. 

We  give  the  results  of  several  experiments.  In  all  the 
examples  we  take  the 

initial  guess  for  aQ(x,y>  =  0.5625. 

In  the  following  we  give  results  in  terms  of  the  relative  L2 
errors  for  a0  at  a  particular  iteration  that  are  calculated  as 
fol lows 

2  <fo  <*0calc  -  »0tst)2  dx  dy  )l/2 

Relative  L  error  ■  — “ - - 

»*0tst  V2(S2) 


i. 


*0tst(x,y)  *  1.0  +  0. 25*cos(2nx)*cos(2ny) 


iter at  ion 

Relative  L*  error 

0 

0.462 

1 

0.299 

5 

0.143 

10 

0.0894 

15 

0.0684 

_ 35 _ 

_ 0.0428 _ 

Example  2. 


*2+  J2 

a0tst<x'y>  =  1  +  ,V 


iter at  ion 

Relative  L*  error 

0 

0.610 

1 

0.410 

5 

0.213 

10 

0.141 

15 

0.0987 

_ 35 _ 

0.0399 _ 

Example  3. 
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aOtst* 


x  |  0.5  and  y 
otherwise 


|  0.5 


iter at  ion 

Relative  error 

0 

0.350 

l 

0.363 

3 

0.  180 

10 

0.112 

15 

0.0692 

_ 25 _ 

_ 0*3323 _ 

If  data  arc  aval  able  as  a  time  series,  we  wish  to  use  data  at 
each  time  step  to  estimate  the  coefficient  "a”  as  the  system  is 
evolving  in  time.  That  is,  at  each  time  step  new  data  is  available 
to  use  to  estimate  a.  Thus,  as  the  von  Karman  system  is  being 
solved  iteratively,  we  seek  to  estimate  a  at  the  current  step  with 
the  data  using  the  new  information.  Similar  methods  are  used  in 
meteorology  to  incorporate  measurements  into  the  solutions  calcu¬ 
lated  by  meteorological  models  and  are  refered  to  as  forward 
assimilation  techniques.  We  apply  these  methods  to  estimation 
problems  for  the  case  of  the  von  Karman  equations. 

Specifically  we  consider  the  time  approximation  with  a  step 
of  lengh  h 

utt  <uk+l  -  2  uk  +  uk-l>/h2 
and  introduce  an  error  function  w  *  w(a,u)  such  that 


(2. 5) ( i )  L  w  -  uk+1  -  2  uk  +  uk_t  + 

+  h2<A  (a  a  uk>  +  c  CB(uk,uk>,uk3  -  fk>in  Q 
with  initial  conditions 


C2.5Mii> 

and  boundary  conditions 


U (0)  ■  Uq 

uj  *  h  ut+  uq 


uk  ■  »  0  on  9Q  x  (0,  T) 

where  B(u,v>  «  H2<£2)  is  the  solution  of  the  boundary  value 
problem 


a2  B(u,v)  «  Cu,vD. 

We  give  two  approaches  to  the  problem  that  we  consider, 
first  is  to  replace  the  term 


The 


-10- 


*  <  a  ^  uk>  by  a  <  a  a  uk+1>, 
set  w  ■  0t  and  consider  the  elliptic  problem 

uk+i+  (a  a  uk+1  )  ■  2  uk  -  uk-1-  c  h2(CB<uk,uk),uk3  -  f) 

for  the  estimation  of  the  coefficient  a.  This  approach  would  take 
advantage  of  the  existence  and  uniqueness  theory  of  local  solu¬ 
tions.  A  second  formulation  is  weaker  in  the  sense  that  we  chose, 

for  example,  the  operator  L  such  that  L  «  ai2  with  the  error  func¬ 
tion  w  satisfying  the  boundary  conditions 

w  -  §f  -  0. 

In  either  case  given  data  zk  at  time  t  ■  kh  with  k  I  1,  we 
look  for  an  out put -least -squares  solution  to  the  estimation  prob¬ 
lem.  In  the  first  case  we  consider  the  functional 

J<a>  -  lluk+1<a>  -  zk+1ll?  +  P  lt«l|22(n) 

where  the  space  Z  is  the  observation  space,  for  example,  L2(Q). 

In  the  second  case  we  consider  the  functional 

J(a, u; K)  -  llu  -  zk  +  1  II?  +  P  +  K  II*  w<a, u)  ll22 

where  K  is  taken  as  a  penalty  to  force  the  error  w(a,u>  *  0  at 
time  t  =  (k+l)h.  Thus,  a  is  estimated  at  each  time  step.  It  would 
seem  natural  to  initialize  the  minimization  iteration  at  the 
current  time  step  using  the  estimated  coefficient  from  the 
previous  time  step.  As  information  becomes  available,  we  then 
refine  the  estimate  of  a  by  minimizing  the  fit-to-data  functional. 

By  utilizing  the  information  in  this  way  we  may  obtain 
estimation  algorithms  that  are  faster  than  those  using  the  entire 
history  for  matching.  Moreover,  this  approach  lends  itself 
naturally  to  adaptive  estimation  of  parameters.  Such  adaptive 
estimation  would  be  appropriate  for  on-line  identification 
necessary  in  the  implementation  of  feed-back  control  strategies. 

3.  Narrow  plate  models. 

In  work  with  Professor  D.  L.  Russell  C5,&3,  we  have  been 
investigating  models  of  structural  elements  that  are  of  an 
intermediate  character.  Specifically,  not  only  do  they  have  a 
small  thickness  but  their  width  is  small  when  compared  to  their 
length  as  well.  Such  elements  arise  frequently  in  aircraft  wings. 


propellers,  and  fans  blades,  for  example.  We  give  the  essentials 
of  the  formulation. 

Let  a  region  ft 

ft  -  < (x, y, z )  :  0  |  x  $  L,  kt(x)  f  y  i  k2(x), 

-  h<x,y)  j  z  J  i  h<x,y)  > 

be  occupied  by  an  elastic  isotropic  body  satisfying  the  small 
deformation  gradient  assumption.  We  assume  that  the  body  is 
clamped  along  the  line  x-z  =  0,  kj(x)  i  y  f  k2<x).  We  also  take 
the  functions  k  j  C  •  >  and  k2(  •)  to  be  of  the  form 

k  ^(x)  *  p  x  -  k(x)  and  k2(x)  =  p  x  +  k<x) 
where  k  is  a  piecewise  continuous  function.  For  our  model  we 
specify  the  displacement  relations  as  a  particular  form  of  the 
Mindl in-Timoshenko  relations  where  u,  v,  and  w  are  the  displace¬ 
ments  in  the  x,  y,  and  z  directions,  respectively. 

u(x,y,z,t)  =  z  £>(x,t>  +  (y-px)  <p(x,t)3 

v<x,y,z,t>  =  z  C*(x,t)  +  (y-px)  y(  x ,  t )  3 

w(x,y,z,t)  =  Z  <x,  t )  +  (y-px)  e(x,t)  +  (y-px)2  o(x,t) 

Simpler  models  that  may  be  useful  are  obtained  by  omitting  the 
quadratic  term  in  (y-px)  in  the  expression  for  w  and  the  linear 
terms  in  u  and  v.  These  displacement  assumptions  give  rise  to  the 
system  in  which  h  and  k  are  constants  and  D  is  the  flexural 
rigidity  and  6  is  the  twisting  modulus. 

§p  *  +  Gkh  (*  +  Z  -  pq)  -  (2kD  ^  I)  M1 

C3.n  ^  +  Qkh  <*  +  "  Is-  <  kD  I ST  *>  “  f2 

2  Pkh  2^2  Z  -  §g.  (Qkh  <  ¥  +  §?  Z  -p«>)  ■  f3 

Ph 1  y  |^2  a  +  Qkh(  *  +  €  -  p(l  +  Isj-  Z  -  P  «))-  f^.  (^5-Iy  §g-  «)  «  f4 

These  equations  are  accompanied  with  clamped  boundary  conditions 
at  x  *  0  and  free  conditions  at  x  =  L.  Clearly,  by  including 
additional  terms,  we  can  accomodate  more  general  motions. 

Further,  by  considering  a  family  of  coupled  systems  of  the  form 
(3.1),  it  is  possible  to  model  corrigated  and  fan-like  struc¬ 
tures. 

There  are  several  advantages  in  considering  the  single  space 
dimensional  model  (3.1)  in  place  of  two  dimensional  plate  models. 
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Th»  first  is  the  dramatic  computational  speedup  and  savings  in 
storage  that  result  from  replacing  a  2  dimensional  problem  with  a 
1  dimensional  system.  Secondly,  the  narrow  plate  model  incorpor¬ 
ates  the  boundary  behavior  into  the  coefficients  of  the  equations. 
In  so  doing  the  regularity  of  the  solutions  may  be  deduced  much 
easier  than  in  the  2— D  model.  The  2— D  model  regularity  is  deter¬ 
mined  by  smoothness  of  the  boundary  where  in  the  1 -D  model  the 
coefficients  are  obtained  by  integrating  the  functions  that  define 
the  boundary. 

This  regularity  allows  us  to  obtain  a  very  nice  convergence 
theory  for  finite  element  approx imat ions,  in  fact  obtaining  rate 
of  convergence.  The  regularity  theory  for  the  2-D  model  is  much 
more  complicated.  Mhile  we  can  obtain  convergence,  a  rate  is 
difficult  to  obtain  except  for  smooth  boundaries.  Such  results  are 
important  since  we  plan  to  study  the  relation  between  the  cross- 
sectional  shape  of  narrow  plate  and  its  the  natural  frequen¬ 

cies  and  mode  shapas.  Incidently,  results  indicate  that  torsional 
frequency  information  is  critical  for  determining  the  shape  of  the 
plate  so  that  it  is  important  to  have  a  theory  that  accounts  for 
more  than  beam-type  bending. 

Me  are  also  considering  the  use  of  spectral  data  for  the 
identification  of  material  parameters  in  narrow  plate  models  C53. 
Accordingly,  we  have  developed  a  theory  using  a  weighted  least 
squares  fit-to-data  criterion  that  supports  algorithms  for  the 
estimation  of  such  parameters.  Qur  theory  applies  whether  or  not 
the  multiplicity  of  natural  frequencies  is  known.  Moreover,  we 
are  determining  the  approximat ion  properties  and  stability  with 
respect  to  perturbat ions  in  the  spectral  data.  Further,  we  are 
testing  our  algorithm  with  laboratory  data  obtained  from  a 
particular  sample  for  which  the  assumptions  of  the 
narrow  plate  model  hold.  The  goal  is  to  identify  the  shape  of  the 
sample.  Me  have  obtained  encouraging  initial  results  that  are 
consistent  with  expected  qualitative  results  predicted  by  the 
implicit  function  theorem.  Investigations  in  this  direction  are 
cont inuing. 


4. 


A  system  of  connected  beams 


During  this  project  we  have  been  conducting  an  investigation 

concerning  the  estimation  of  coefficients  or  the  design  of  beam 
systems  C7D.  In  particular,  we  have  recently  been  studying  the 
following  simple  problem.  Suppose  two  beams,  1  and  2,  of  length 
and  ?2 1  respectively,  are  Joined  at  an  angle  oc  to  one  another. 
The  other  end  of  beam  1  is  clamped  and  the  remaining  end  of  beam 
is  free.  In  response  to  a  force  that  is  perpendicular  to  the 
plane  containing  the  beams  are  deformations  wj  and  w2  of  beams  1 
and  2,  respectively,  and  a  rotation  6  of  twist  along  beam  1.  The 
key  point  for  the  model  we  have  been  investigating  is  that  the 
angle  between  the  beams  remains  fixed.  For  the  dynamic  problem 
with  oc  =  90°,  we  study  the  following  system  of  hyperbolic 
equations  with  coupled  boundary  conditions  in  which  x  and  y 
represent  local  coordinates  for  beams  1  and  2  respectively. 

<a  wxx>xx  “  fl  in  <°»  »1> 

C4. 1 )  C  i  )  -<b  ex)x  *  g  in  <0,  ) 

<c  w2yy * yy  ~~  ^2  ^  f?2) 

with  boundary  conditions 

Wj  <0)  ■  Wjx  CO)  *  0 
SCO)  *  0 

at  the  clamped  end, 

C4. 1) Cii)  wlxx( *!>  =  0 

C4.  1)  Ciii)  Wj  C  Pj)  -  w2C0) 

C4.  1)  Civ)  9C  i?t)  *  w2y CO) 

C4.  1 )  C v)  Cb  ©X)CI?1)  -  <c  w2y y )  C 0 ) 

C4.  1 )  C  vi  )  Ca  wlxx)x(Pj)  -  Cc  w2yy)yC0) 

at  the  Junction,  and 

w2yy<p2)  "  <c  w2y y ) y (  92 >  =  0 
at  the  free  end.  Equations  C4.1)  are  associated  with  a 
potential  energy  functional 

C4.2)  PCW},  9,  w2)  ■  CaCx)  w^xxCx)  +  bCx)  0^Cx)D  dx  + 

+  Xq2  w|yyCy)  dy  -  2  CftCx)  w^x)  +  gCx)  0Cx)D  dx  - 

”  2  J 0  f2<y)  w2Cy)  dy. 

A  unique  solution  of  the  initial  boundary  value  problem  C4.1) 
may  be  shown  to  exist  using,  for  example,  standard  Galerkin 
arguments  with  the  Hilbert  space 
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V  -  <v  »  <wlf©tw2)  m  H2<0,  Jj)  x  Hl(Of  >  x  H2<0,  02>  s 
wt<0)  ■  wlx<0>  •  Of  0<O>  *  0, 

$<!>!>  *  w2y<0)f  WjCJj)  »  w2<0)>. 

To  solve  (4.1)  numerically  with  a  comforming  method,  the 
space  of  basis  functions  should  satisfy  the  essential  boundary 
conditions  given  in  the  definition  of  the  space  V.  However,  it  is 
difficult  to  enforce  the  essential  boundary  conditions  at  the 
Junction  on  the  basis  elements  that  are  used  to  construct  the 
finite  dimensional  approximating  space,  as  one  does,  for  example, 
with  the  clamped  boundary  conditions.  We  are  motivated  by  a 
desire  to  develop  and  study  methods  that  can  be  generalized  to 
more  complicated  structures.  Hence,  our  approach  is  to  use  basis 
functions  that  are  nonconforming  at  the  Junction  in  the  sense 
that  they  do  not  satisfy  the  essential  boundary  conditions  at  the 
Junction.  For  a  static  problem 

these  conditions  are  enforced  as  constraints  on  the 

discrete  version  of  the  problem  to  minimize  the  potential  energy 

functional  and  a  penalty  method  is  used  to  impose  the  conditions. 

For  the  dynamic  case  we  consider  we  consider  an  initial  boundary 
value  problem  having  a  potential  energy  functional  that  is 
obtained  through  a  penalization  of  the  functional  (4.2). 

To  give  a  penalization  formulation  we  define  by  W  the  space 

W  ■  <v  ■  <wlt0,w2)  m  H2<0,  x  H1(0,  (?i>  x  H2<0,  l>2>  : 
wj  (0)  *  Wjx<0)  *  0,  0(0)  =*  0>. 

and  the  mapping  G:  V  R2  by 


G(w)  =  fwj  <  -  w2 CO)  1 

L«P<  -  w2y<0)  J. 

In  order  to  obtain  convergence  results,  apparently  it  is  useful  to 
define  approximat ing  problems  in  terms  of  a  penalized 
potential  energy  functional  to  incorporate  the  constraint  at  the 
Junction.  For  example,  one  such  functional  is  given  by 


J^Cu)  =  PCu)  +  K  |G(u)  |2 


where  «  e  <0,1).  We  point  out  that  apparently  regular izat ion  is 
not  necessary  in  this  case  as  in  the  static  case  C73. 

We  are  currently  considering  estimation  problems  for  this 
structure  based  on  the  potential  energy  functional  (  •  ' .  In  this 
work  we  formulate  a  sequence  of  estimation  problems  where  in  each 
case  the  underlying  system  is  one  obtained  by  the  penalizing  the 
condition  at  the  Junction  thereby  absorbing  the  Junction  condition 
into  a  regularized  functional.  Hence,  a  sequence  of  estimation 
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problems  is  defined  on  the  regularized  approximating  systems 
obtained  from  the  penalty  approximating  scheme. 
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